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CO ■ An explicit parameterization for the state space of an n-level density matrix is given. 

^^D ' The parameterization is based on the canonical coset decomposition of unitary matrices. 

^p , We also compute, explicitly, the Bures metric tensor over the state space of two- and 

JN ■ three-level quantum systems. 
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In recent years, the Riemannian Bures metric [1] has become an interesting subject for the 



1 Introduction 

Sh| understanding of the geometry of quantum state space. It is the quantum analog of Fisher 

information in classical statistics, i.e. in the subspace of diagonal matrices it induces the 
statistical distance [2] . The Bures measure is monotone in the sense that it does not increase 
^^ . under the action of completely positive, trace preserving maps [3]. It is, indeed, minimal 

5-H ' among all monotone metrics and its extension to pure state is exactly the Fubini-Study metric 

[3] . The Bures distance between any two mixed states pi and p2 is a function of their fidelity 
Fipi,P2) [4,5] 



dBipi,P2)^^j2-2^F{pi,p2), F{pi,p2) 



Tr a/ VPiP2\/pi" 



(1) 



Fidelity allows one to characterize the closeness of the pair of mixed states pi and p2, so, it 
is an important concept in quantum mechanics, quantum optics and quantum information 
theory. An explicit formula for the infinitesimal Bures distance between p and p + dp was 
found by Hiibncr [6] 

dB(p,P + dp) =-)^^-_^^^, (2) 
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where Xj and |j), (j = 1, 2, • • • , n) represent eigenvalues and eigenvectors of p, respectively. 
Dittmann has derived several explicit formulas, that do not require any diagonalization pro- 
cedure, for Bure metric on the manifold of finitc-dimcnsional nonsingular density matrices 
[7, 8], for instance 



dB{p,p + dp) = -Tr 



dpdp + — (dp - pdp){dp - pdp) 



(3) 



and 

dB{p,p + dp) =jTt dpdp + (i_T^p3) (dp - pdp) (dp - pdp) 

+ (i-Tr'p-') (^P ~ P~^dp)(dp - p-^dp) , 

for nonsingular 2x2 and 3x3 density matrices, respectively. 

The probability measure induced by the Bures metric in the space of mixed quantum states 
has been defined by Hall [9] . The question of how many entangled or separable states there are 
in the set of all quantum ststes is considered by Zyczkowski et al in [10, 11]. Sommers et al [12] 
have computed the volume of the (n^ — l)-dimensional convex set and (n^ — 2)-dimensional 
hyperarea of the density matrices of an n-level quantum system. In a considerable work, 
Slater investigated the use of the volume elements of the Bures metric as a natural measure 
over the (n^ — l)-dimcnsional convex set of n-level density matrices, to determine or estimate 
the volume of separable states of the pairs of qubit-qubit [13, 14] and qubit-qutrit [15, 16]. 
Very recently [17] , Slater made use of the Bloore parameterization [18] of density matrices in 
order to obtain the Hilbert- Schmidt volumes of separable subsets for the two qubit system. 

The state space of an n— level quantum system is identified with the set of all n x n Hermi- 
tian positive semidefinite complex matrices of trace unity, and comprise (n^ — l)-dimensional 
convex set. Due to considerable interest in the use of density matrices, a lot of work has 
been devoted to describe and parameterize density matrices. Any density matrix of an n-level 
system can be expanded in terms of orthogonal generators Aj of SU{n) as [19] 

1 / n(n-l)_, ->\ 

P=-(/„ + ^^.-AJ (5) 

where ly — (i^i, 1^2 ■■ ■ , f^n^-i) is a real vector and Lie algebra generators Xi are normalized 
as Tr(AiAj) = 2Sij. The above representation is the generalization of the Bloch or coherence 
vector representation for two-level systems and gives one of the possible descriptions of a state 
on the basis of the actual measurements which is an important task both from experimental 
and theoretical viewpoint [20]. The region of Bloch vector v which represents a physical 
density matrix have been found in [20, 21]. An investigation of the geometrical aspects of the 
Bloch vector space from the spherical coordinate point of view is also made by Kimura et al 
in [22]. 

Boya et al [23] have shown that the mixed state density matrices for n-level systems can 
be parameterized in terms of squared components of an (n — l)-sphere and unitary matrices. 
By using the Eulcr angle parameterization of SU{3) group [24], 

U = g«-*'3ag«>'2/3giA37giA5egiA3agiA2&giA3CgiA80/V3 /g\ 

where Xi are Gell-Mann matrices, Byrd et al [25] have presented a parameterization for eight- 
dimensional state space of three-level system as 



P 



^ VDV\ (7) 
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where Z? is a diagonal density matrix (with two independent eigenvalues), and V G SU{3) is 
given by 

They also gave Bures measure on the space as the product of the measure on the space of 
eigenvalues and the truncated Haar measure on the space of unitary matrices. An Euler 
angle-based parameterization for density matrices of four-level system is also introduced in 
[26]. A generalized Euler angle paramctrization for SU{n) and U{n) groups has given by 
Tilma et al [27, 28]. Tilma et al [29] have also used the parameterization for four-level system 
(two qubit system) in order to study the entanglement properties of the system. 

In a comprehensive analysis [30], Zyczkowski et al analyzed the geometrical properties of 
the set of mixed quantum states for an arbitrary n-level system and classified the space of 
density matrices. Di^a [31] has provided a parameterization for general Hermitian operators 
of n-level quantum systems. The parameterization is based on the factorization of n x n 
unitary matrices [32] and may be used either for Hamiltonian or density matrices. In [33], 
the authors have shown that the space of two qubit density matrices (four- level systems), can 
be characterize with 12-diincnsional (as real manifold) space of complex orthogonal group 
SO{4:,C) together with four positive Wootters's numbers [34], where of course, the normal- 
ization condition reduces the number of parameters to 15. 

By using the definitions (4) and (7), Slater has computed the Bures metric for the eight- 
dimensional state space of three-level quantum systems [35] . He showed that all entries of the 
8x8 matrix tensor are independent of the Euler angle a, and the matrix tensor decomposes 
into a 6 X 6 block and a 2 x 2 one, in corresponding to the six Euler angles of unitary matrix 
V and the two independent eigenvalues of diagonal matrix D. 

In this paper we consider a canonical coset parameterization for density matrices of an 
n-lcvel quantum system. The parameterization is based on the coset space decomposition of 
unitary matrices [31]. This parameterization, as well as the Euler angle parameterization do, 
eliminates any over-parameterization of the density matrix. It also provides a factorization of 
the Bures measure on the space of density matrices as the product of the measure on the space 
of eigenvalues and the truncated Haar measure on the space of unitary matrices. We give 
explicitly, the parameterization for two- and three-level density matrices, and by application of 
Dittmann's formulas, the Bures metric over the spaces of these quantum systems are explicitly 
computed. It is shown that the coset parameterization gives a compact expression for the 
all metric elements. The analytical simple expression obtained for the Bures metric of the 
three-level system enable us to use this parameterization for the problem of Bures metric over 
the space of two qubit system. 

The paper is organized as follows: In section 2, the coset space parameterization of an 
n-level density matrix is introduced. We give also, explicitly, the parameterization of the 
density matrices of the two- and three-level systems in this section. In section 3 we compute 
explicitly the Bures metric of the two- and the three-level systems. The paper is concluded 
in section 4 with a brief conclusion. 

2 Canonical coset parameterization of density matrices 

In this section we review some properties of the set of density matrices of an n-level quantum 
system and, by using a canonical coset parameterization for nxn unitary matrices, we present 
a coset parameterization for n-level density matrices. The state space of an n— level quantum 
system is identified with the set of all n x n Hermitian positive semidefinite complex matrices 
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of trace unity, and comprise (n^ — l)-dimensional convex set Mn- The total number of 
independent variables needed to parameterize a density matrix p is equal to n^ — 1, provided 
no degeneracy occurs. 

Let us denote the set of all diagonal density matrices of an n- level system with I?„. An 
arbitrary element D G I?„ can be written as 

n 

i? = diag{Ai,A2,--- ,A„}, 0<A, <1, ^A, = 1. (9) 

This means that the set of all diagonal density matrices forms an (n — l)-dimensional simplex 
Sn-i . A generic density matrix in an arbitrary basis can be obtained as the orbit of points 
D G Vn under the action of the unitary group U{n) as 

p = UDU\ (10) 

Let iJ be a maximum stability subgroup, i.e. a subgroup of U{n) that consists of all the 
group elements h that will leave the diagonal state D invariant, 

hDh'^ = D, he H, DeVn, (11) 

that is, H contains all elements of U{n) that commute with D. For every element U G U{n), 
there is a unique decomposition of U into a product of two group elements, one in H and the 
other in the quotient G/H [36], i.e. 

U^nh, UeU{n), heH, VteU{n)/H. (12) 

The above decomposition implies that the action of an arbitrary group element U G U{n) on 
the point D G I?„ is given by 



P 



= UDU^ = nhDh'^9) = ^DVtK (13) 



This means that in order to characterize the space Al„, it is sufficient to consider the or- 
bit of points D G Vn under the action of the quotient U{n)/H. Since I?„ consists points 
with different degree of degeneracy, the maximum stability subgroup will differ for differ- 
ent D £ Vn [30]. Let nii denotes degree of degeneracy of eigenvalue \i of matrix D. 
This kind of the spectrum follows that D remains invariant under the action of arbitrary 
unitary transformation performed in each of the ?7ii-dimensional eigensubspace. Therefore 
H = U{mi) ® U(m2) ® • • • U{mk) is maximum stability subgroup for D, and the quotient 
space U{n)/H is a complex flag manifold 

U{n) , , 

•^ = 777 ^ TTT ^ 777 T , mi + m2 + ■ ■ ■ + mk — n. (14) 

C/(mi) (g) [/(ma) (g) • • • L/(TOfc) ^ ^ 

Two special kinds for the degeneracy of the spectrum of D are as follows: i) Let D represents 
the maximally mixed state p* = diag{i, —,■■■,—}. In this case the stability subgroup H 
is U{n), and the orbit of point p^ is only one point, i.e. p ^ p^. ii) On the other hand 
if the spectrum of D is non-degenerate, then the stability subgroup is n-dimensional torus 
T" = [/(l)"^", and the orbit of the point D is 

p = nDn\ n g u{n)/T". (15) 
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The maximal torus T" is itself a subgroup of all maximum stability subgroups, therefore the 
orbit of points D e I?„ under the action of quotient U{n)/T^ generates all points of the space 
AAn- The diagonal matrix D is defined up to a permutation of its entries and, one can divide 
the simplex 5„_i into n\ identical simplexes and take any of them. Each part identify points 
of iS„_i which have the same coordinates, but with different ordering, and can be considered 
as the homomorphic image of simplex 5„_i relative to the discreet permutation group P„, i.e. 
Sn-il Pn- Therefore the points of Al„ can be characterize as the orbit of diagonal matrices 
D e Sn-i/Pn under the action of quotient U{n)/T". 

Further insight into the space of density matrices can be obtained by writing to the 
arbitrary element U € U{n) as [36] 

u^ni:'^n^^\---ni"^n^;'\ (i6) 

where 17^"^ e T", and 

(„) U{k) T-'- fc-2... n (171 

Comparing this with the decomposition U = Clh, where h G T", leads to the following 
decomposition for an arbitrary clement fl of quotient U{n)/T" 

n^ni-'>n^S---ni"\ (is) 



A typical coset representative ilj^ can be written as 

(„) _ / SUik)/U{k-l) 
^h " I qT 



^ ' (19) 



J-71—k 



where O, O^ and In-k represent, respectively, the k x {n — k) zero matrix, its transpose and 
the (n — fc) X (n — k) identity matrix. The 2(fc — l)-dimensional coset space SU{k)/U{k — 1) 
have the following k x k matrix representation [36] 



SU{k)/U{k-l) 



;/BW[SW]t 



/rR(fc)itB(*:) I- ' 



VlB 



V[^"°'l^^'^'' I (20) 



cos 



V[BW]tB(fe) 



where B^'^' represents (fc — 1) x 1 complex matrix and [i?''^-']^ is its adjoint. In the following 
we consider the n — 2,3 cases, explicitly 

2-1 Two-level system 

We begin by giving the coset parameterization for a two-level quantum system. Let us consider 
a diagonal two-level density matrix D = diagjAi, A2}. Every h E H = T^ leaves the density 
matrix D invariant i.e. HDh^ — D for h E H. Any group clement U G U{2) can be 
decomposed, uniquely, a,s U — flh where ft E C/(2)/T^ and h E T^ [36]. The coset space 
with respect to the stability subgroup H = T^ will provides the unitary transformations to 
construct a generic density matrix p as the orbit of diagonal matrix D, i.e. 

p = nDn\ nEU{2)/T^. (21) 
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A typical coset representative in the coset space U{2)/T'^ is 

n^ ( ^^^'^ e^-^sina \ ,^^, 

y — e~*'^sina cos a y ' 

where a,(f) are real. The range of parameters Ai, A2 can be determined as follows. The set 
of all diagonal 2x2 matrices D forms a 1-dimensional simplex Si , which can be divided into 
two identical parts 0<Ai<^,^<Ai<l (see figure la). It can be easily seen that each 
part can be obtained as the orbit of the other part under the action of the group element 
il ~ n{a ~ (p ~ ^) which is an element of the coset space U(2)/T'^. This means that we can 
easily consider the diagonal matrix D as 

D=(""t' .',V 0<e<l (23) 

\^ sm^ e J ' ~ ~ 4 ^ ^ 

With this parameterization any 2x2 density matrix can be written explicitly as 

/ sin^ q; sin^ + cos^ a cos^ 6* —ie*"^ sin 2a cos 20 \ , , 

^^\ -ie-'"^ sin 2a cos 26* sin^ a cos^ 6* + cos^ a sin^ e* j ' *- '' 

2-2 Three-level system 

The density matrices for the three-level system comprise eight-dimensional convex set. Let 
D = diagjAi, A2, A3} be a 3 x 3 diagonal density matrix of a three-level system. The density 
matrix D is invariant under the action of every group element h E H = T^, i.e. HDh^ = D 
for h E H. The coset decomposition U ~ ilk where il E U{3)/T^, h E T^ provides a 
parameterization for a generic p as 

p = UDU^ = nhDh^n^ = r^L*!]!'. (25) 

On the other hand the decomposition 

u^nf^nf'>n^^\ (26) 



with 



follows that 



nf'> e C/(3)/C/(2) ® C/(l), 

rif e (C/(2) ® t/(l))/(t/(l) ® C/(l) C/(l)), (27) 



n^i^ ec/(i)®c/(i)®c/(i), 



n^n^^''n^^K (28) 



(3) (3) 

The coset representatives ilj and fl^ can be parameterized respectively as 

■,(3) 




n^2 ^ ^ \ -e-^Psina cosa , (29) 



and 



l-h|i(cos/3- 1) ^e*('^i-'''=^)(cos/3-l) ^e*'^! sin/3 

^3'^=! ^e-*('^i-'^^)(cos/3-l) l + §|(cos/3-l) ^e'^^sin/3 I, (30) 



^e-'^^ sin 13 -| 



e-''^! sin/3 -^e"''''^ sin/3 cos/3 
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where /3 = ^/|3f + f3^ . 

The two-dimensional simplex 52 of eigenvalues of p is divided into 3! parts (see figure 
lb). Let us take one of the parts (e.g. the shaded one) for illustration. It is easy to see 
that all other parts can be obtained from this one by applying the permutation group P3. 
On the other hand the elements of the permutation group P3 can be obtained from the coset 
representative (28) up to a phase as 






\(ii = f ,/32 = 0,7/^1 = f )f^r (a = 4>=i)= ^(123), 
\Pi = 0,/32 = f ,^2 = %)^f{a = 0= f ) = z(321). 



/32- 


QW^'\a 


- 0) = [Id), 






/32- 


Q)9.^^\a 


-^=f)- 


2(12), 




f,/32 


= 0,1/^1 = 


f)"f (a = 


= 0) = 


z(13), 


0,/32 


= f>2 = 


f)«f (a = 


= 0) = 


2(23), 



^i^)| 



Therefore the ranges of the eigenvalues of D can be determined as -^ < Ai < 1, < A2 < 5 
and < A3 < i, or cquivalcntly, one can parameterize the diagonal matrix D as 



C0S2 01 











sin^ 01 cos^ 62 











sin^ 9i SVC? 02 



D= \ sin^ 6*1 cos2 6*2 , (31) 



where < 6*1 < cos"! -i^, f < ^2 < f • 

3 Bures metric 

In this section we calculate the Bures metric of the two- and the three-level quantum systems. 
We will use the canonical coset parameterization of the density matrices introduced in the 
last section. 



3-1 Two-level system 

In the two-level systems by using Eqs. (3) and (21) we get 

dB{p,p + dp) = iTr [{dDf + ^(X - DfidDf 
+ E.<, [%(("^df7),,(f^tdf7),,)], 



where 2 is unit matrix and 

1 



bi' 



\D\ 



[{Du - D,,f{Du + D,, - DuD,, - \D\ - 1)] . 



(32) 



(33) 



By defining the ordering {9, a, (j)} for coordinates, the corresponding Bures metric tensor takes 
the following diagonal form 



1 

cos2 29 














isin^2Q:cos2 26' 



(34) 
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(0,1) 



(a) 



(1,0) 



(0,1,0) 




1,0,0) 



ro i i 

^", 2 ' 2 



(5,0,i) 



Figure 1: (a) Onc-dimcnsional simplex Si of diagonal density matrices of two- level systems, 
(b) the two-dimensional simplex 52 of diagonal density matrices of three-level systems. The 
simplex 5„ can be decomposed into n\ parts. The parts can be transformed to each other by 
applying the elements of the permutation group P„. 



3-2 Three-level system 

In the three-level system, by using (25) in (4) we have 



dB{p,p + dp) = |Tr [{dDf + (^_ jri^3) 



{{l-Df{dDf + \D\{I - D-^f{ADf 



where 



^ij 2(1 — TrD-^) W^ii jj) \^ii ' -^jj i^ jj -Lri/ zj 

+ \D\{D-^ - D-^nn-:' + Djl - Dr^D-^ 1)] . 



(35) 



(36) 
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Now defining the ordering of coordinates as {^i, ^2, a, 07/9ij/52, V'l: V'2}, the corresponding 
symmetric matrix for metric tensor takes the fohowing form 



/ 1 

















\ 


sin^ 9i 






















9aa 





5a/3i 


9af32 


9atpi 


5ai/'2 






9<P<P 


9<pl3l 


94>p2 


94>4>i 


9M2 








5/3i/3i 


5/3i/32 

5/32/32 


5/3iVi 

902->Pl 

9i'ii'i 


5/3iV2 

5/32 V'2 

9il> ii>2 


\ 












9i>2lp2 / 



(37) 



The 8x8 matrix tensor decomposes into a 2 x 2 block and a 6 x 6 one, in corresponding 
to the two independent eigenvalues ^1,^2 and six coset parameters a, 0, /3i,/32, V'l: V'2- After 
some analytical calculations, we can obtain the following expression for the matrix elements 
of the 6x6 block 

• Qnrv = —Tl2 



• 9a4, = 

• 5a/3i 

• 5q/32 

• 5a^l 

• 9aif>2 

• 500 = 







2Ti2/32COS7^^-g- 

_/3i 



/32 



5a/3i 



2Ti2/3i/32(72sin7(^y 



U2 9aipi 

^Ti2 sin^ 2a 



50/32 
50-01 
50-02 

5/3i/3i 



j-Ti2/32sin4Q: 1^ ^ 



■fc50/3i 



\Ti2(3i sin 2a [/3iW2sin2Q! + 2/32C/2 cos 2q: cos 7] 



= ^ri2/52sin2Q: [/32W1 sin2Q! — 2/3iC/i cos2acos7] I — g^ 



= -4Ti2/3| [1 - sin^ 2a sin^ 7] , ^ 

— Tia X^ sin^ a + y,^ cos^ a — XVi sin 2a cos 7 
-r23 [X^ cos^ a + V^ sin^ a + XFi sin 2a cos 7 



5/3i/32 = 4Ti2/3i/32 [1 - sin^ 2a sin^ 7] 



-Ti3 

-T23 



X(yiCos2a + y2sin^a) - 5(ViV2 + X^) sin 2a cos 7 
X(yisin^a + y2C0s2a) + i(ViV2 + X^) sin 2a cos 7 
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ff/3ii/>i = ^Ti2l3il32 [2P2U2 sin^ 2a sin 27 - /3iT4^2 sin 4a sin 7] 
+ 5(Ti3 - T23)Pi sin2asin7 [C/2X + ViY] (siE^) 

9M2 = -Ti2(32 [2/3iC/isin2 2asin27 + /32VKisin4asin7] (^) 
■5(^13 - T23)/32 sin2asin7 [U^Vi + XF] (^' 



5/32/32 



-4^12/3? 1 - sin^ 2a sin^ 7 



/3 



-^13 
-2I23 



X^ cos^ a + V2 sin a — XV2 sin 2a cos 7 
X'^ sin^ a + Vj^ cos^ a + XV2 sin 2a cos 7 



5/32^1 = TuPf [2P2U2 sin2 2a sin 27 - (3iW2 sin 4a sin 7] (^) 



1(^13 - r23)/3i sin2asin7 [[72^2 + ^>^] (^ 



5/32^2 = Ti2(3i(32 [2/3iC/i sin^ 2a sin 27 + /32W^i sin 4a sin 7] 



-i(ri3 - T23)/32 sin2asin7 [(7iX + V2Y] (^) 
ffVi^i = -T12PI [4/3|C/|(l - sin^ 2acos2 7) 

+/3^ Wa^ sin^ 2a + 2/3i/32 U2W2 sin 4a cos 7 



-TiiPl [Ui cos2 a + r2 sin2 a + U2Y sin 2a cos 7] ( ^ 



sin f3 



-T23.(il [Ul sin^ a + Y"^ cos^ a - L/2>^ sin 2a cos 7] 
5ViV2 = -Ti2(3iP2 [4/3i/32C/i 1/2(1 - sin^ 2acos2 7) 

-(31P2W1W2 sin^ 2a - sin4acos7(/3|J72VFi - PIU1W2)] , -g 
+ri3/3i/32 [5^(?/i sin^ a + U2 cos^ a) + i sin2acos7(C/iC/2 + Y^)] 
+T23/3i/32 [i^(?/i cos2 a + C/2 sin^ «) - ^ sin2acos7(C/iC/2 + Y^)] 



sin f3 
sin /3 



+PiW{ sin^ 2a - 2/3i/32C/i VKi sin 4a cos 7 



■ /3\4 



(^) 



-ri3/3| [t/2 gjjj2 ^ _^ y2 ^Qg2 q, + ^^y gin 2a cos 7] f ^ 
-123/31 [Uf cos2 a + y2 sin2 ^ - UiY sin 2a cos 7] (^ 
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In the above equations we have used the foUowing definitions 



7 = </>- Vi +"02, 



and 



Tl2 

T23 



Y" (cos^ 9i - sin^ 61 cos^ 9; ^ 

^ (C0S2 01 
2 



- sin 



'2) 

,2 " ^2 



and 



-1 
2 



C/i 



ii sin^ 62 
sin^ 6*1 cos2 6*2(1 + 3sin2 6'i 



"o I (1-sin'^ fli cos'' e2)(l+cos'' Bi cos"^ 82) 

"^ ' cos* 6*1 cos* 6(2 (cos^ Si+sin* 9i sin^ 62 cos^ 
o , (1-sin^ 6)1 sin^ 92)(cos^ ^2+5111" 9i cos* fli sin* $2) 
•^ ^" sin^ 6*1 cos* 81 cos* e,fcos2 61 +sin* Si sin^ ^2 cos^ 62) 



l2COS^e2) 



(i— sm Wi sm fc'2j(cos 6'2-f 
sin^ 6*1 cos* 81 cos* 82(0032 61 
rns 81 



COS 81 



sin*^ 81 sin* 82 cos^ 82 



fi + f|cos/3, 



t/2 



V, - '^i 



/32^ 



^ COS l3, 



/sin/3\ T. _ Pi , /3? /sin/3\ 



VKi = (1 + cos/3) + 2|i-(l - cos/3), W2 = {1 + cos/3) + 2||(1 - cos/3), 



X 



_ /3i & 



sin f3 \ 



^ =%%(!- cos /3). 



(38) 



(39) 



(40) 



It should be stress that all elements gij are simply products of two independent functions, one 
of the coset space parameters, and the other of the spherical angles 6i, 62- It is also worth to 
note that the three angles 0, ■01 and ■02 appear only in the form 7 = (/> — V'l + "02 ■ 



4 Conclusion 

We present a coset parameterization for density matrices of an n-level quantum system. The 
parameterization is based on the canonical coset decomposition of unitary matrices. By 
using the parameterization for two- and three-level quantum systems, the Bures metric over 
the state space of these systems are computed explicitly. We show that in the canonical 
coset parameterization the symbolic expression for all tensor elements can be obtained. The 
problem of computing the Bures metric of a two-qubit (four-level) system, which is important 
for calculations involving entanglement, is also under consideration. 
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